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In this work, a new model in kinetic gas theory for deriving the Maxwellian Velocity 
Distribution (MVD) is proposed. We construct an operator that governs the discrete time 
evolution of the velocity distribution. This operator, which conserves the momentum and 
the energy of the ideal gas, has the MVD as a fixed point. Moreover, for any initial out- 
of-equilibrium velocity distribution, it is shown that the gas decays to the equilibrium 
distribution, that is, to the MVD. 
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1. Introduction 

The Boltzmann-Gibbs distribution (BGD) and the Maxwellian velocity distribu- 
tion (MVD) are two fundamental distributions in statistical physics. The BGD 
represents the distribution of the states of a system in equilibrium and the MVD 
represents the distribution of velocities in an ideal gas in equilibrium. Although 
they can be directly derived from the equilibrium statistical theory, nowadays there 
was not an unified framework showing that both distributions can be seen as fixed 
points of two different nonlinear operators acting on the space of distributions. Ev- 
idently, it must be cited the Boltzmann theory [1] as a successful attempt in this 
direction explaining the MVD from a microscopic point of view applied to kinetic 
gas theory. 

A scheme inspired in economic systems that has recently been proposed[5] to ex- 
plain the attractivity, and then the ubiquity, of BGD reads as follows: Let p(m)dm 
denote the PDF (probability density function) of money in a multi-agent economic 
system, i.e. the probability of finding an agent of the ensemble with money be- 
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tween m and m + dm. Consider now the discrete time evolution of an initial money 
distribution po(m) at each time step n under the action of an operator T, which 
represents the average effect on the system of many random binary interactions 
(with number of the order of the system size) between pairs of agents exchanging 
their money. Thus, the system evolves from time n to time n + 1 to asymptotically 
reach the equilibrium wealth distribution p f (m) , i.e. 

lim T n (po(m)) -> p f (x) = BGD = 8e- &x with 5 =< Po >~ x . 

n— >oo 

In this case, pf(m) is the exponential distribution (BGD) with the same average 
wealth < pf >— <5 _1 than the initial one < pq >, due to the local and total money 
conservation [3... The mathematical properties of operator T have been disclosed 
in Ref.[4]. Hence, this framework not only puts in evidence that the BGD is the 
equilibrium distribution if not that in this case the BGD is asymptotically reached 
independently of the initial wealth distribution given to the system, a point of view 
that to date was possibly lacking in the literature. 

In this work, we extend this perspective to another problem of the same sta- 
tistical nature. Our goal is to explain the ubiquity of the MVD in ideal gases [5]. 
In the next section, we construct an operator T in the space of velocity distribu- 
tions in order to explain the decay of any initial velocity distribution to the MVD. 
Then, in section [3] we proceed to prove the dynamical properties of this operator 
T, in particular that the MVD is the fixed point of the system where the dynam- 
ics asymptotically evolves. In section |4j several examples showing this dissipativc 
behavior are depicted. Finally, our conclusions are given. 

2. The new model for ideal gases. Operator T 

Consider an ideal gas with particles of unity mass in the three-dimensional (3-D) 
space. As long as there is not a privileged direction in the equilibrium, we can take 
any direction in the space and to study the discrete time evolution of the velocity 
distribution in that direction. Let us call this axis U . We can complete a Cartesian 
system with two additional orthogonal axis V, W. If p n (u)du represents the prob- 
ability of finding a particle of the gas with velocity component in the direction U 
comprised between u and u + du at time n, then the probability to have at this 
time n a particle with a 3D velocity (u, v, w) will be p n (u)p n (v)p n (w) . The particles 
of the gas collide between them, and after a number of interactions of the order of 
system size, a new velocity distribution is attained at time n + 1. Concerning the 
interaction of particles with the bulk of the gas, we make two simplistic and realistic 
assumptions in order to obtain the probability of having a velocity x in the direction 
U at time n + 1: (1) Only those particles with an energy bigger than x 2 at time n 
can contribute to this velocity x in the direction U, that is, all those particles whose 
velocities (u, v, w) verify u 2 + v 2 + w 2 > x 2 ; (2) The new velocities after collisions 
are equally distributed in their permitted ranges, that is, particles with velocity 
(u,v,w) can generate maximal velocities ±U max — zt^/u 2 + v 2 + w 2 , then the al- 
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lowed range of velocities [—U ma x,U ma x} measures 2\U max \, and the contributing 
probability of these particles to the velocity x will be p n (u)p n (v)p n (w) / (2\U max \) . 
Taking all together we finally get the expression for the evolution operator T. This 
is: 



, f f f I'n ■ " U'n \ ' <l'n \ " ! , i . 

Pn+i(x) =Tp n (x) = / / / - ) '= o ==■ Audv&w 

J J Ju 2 +v 2 + 



Pn{u)Pn{v)Pn(w) 

w 2> x 2 2\Ju 2 + v 2 + w 2 



Let us remark that we have not made any supposition about the type of inter- 
actions or collisions between the particles and, in some way, the equivalent of the 
Boltzmann hypothesis of molecular chaos would be the two simplistic assumptions 
we have stated on the interaction of particles with the bulk of the gas. But now 
a more clear and understandable framework than those usually presented in the 
literature appears on the scene. In fact, the operator T conserves in time the energy 
and the null momentum of the gas. Moreover, for any initial velocity distribution, 
the system tends towards its equilibrium, i.e. towards the MVD. This means that 

lim T n (po(x)) -> Pf{x) = MVD (ID case) . 

n—>oo 

Let us proceed to show all these properties in the next section. 
3. Decay to Maxwellian distribution. Properties of T 

First, in order to set up the adequate mathematical framework, we provide the 
following definitions. 

Definition 3.1. We introduce the space L\ of positive functions (one-dimensional 
velocity distributions) in the real axis, 

L+($l) = {p(x):$l^$l + U{0}, |H|<+oo}, 

with norm | • | defined by 

/+oo 
p(x)dx. 
-co 

Therefore for each p e Lfffi), we have ||p|| > 0. 

Definition 3.2. For each x <E 3?, p <E Lf($i), we define the operator T as follows 

J J Ju 2 +v 2 +w 2 >x 2 2vu z +v z +w z 
Therefore we have Tp(x) > for each ieS, and so > 0. 

Remark 3.1. The following relations hold: 

f + OO 

(■)du ) da; 

2 >x 2 
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/+oo / f f \ z-+oo ^+oo / ^ + Vti 2 +f 2 \ 
/ / (-)dudv dx= / / / (-)dx )dudv, 

-oo \J Ju 2 +v 2 >x 2 J J-oo J-oo \J-^u 2 +v 2 J 

f + ( f f f (-)dudvdw) dx = 

J-oo \J J Ju 2 +v 2 +w 2 >x 2 J 



/-TOO /—too r^oo / /— t-VM "t» t» \ 
/ / / {-)dx dudvdw. 

-co J —oo J —oo \J — \/u 2 -\-v 2 +w 2 I 

Theorem 3.1. For any p G Lf($l), we have \\Tp\\ = \\p\\ 3 . 
Proof. Suppose that p G Lf($i) then 

/+oo 
Tp{x)dx = 
-oo 

' ' ' ' ' ' ' p(u)p(v)p{w) 



- ... ^ dudvdw I da; 

2 +v 2 +w 2 >x 2 2yu + v 2 + w 2 J 

,+oo r +oo r +oo / r + Vu 2 +v 2 + W 2 p ( u ) p ( v }p( w ) \ 



/+oo r+oo r+oo I r- 
-oo J-oo J-oo \J-Vu 2 +v 2 +w 2 2\Ju 2 + V 2 + W 2 

/+oo r+oo r+oo 
/ / p(u)p(v)p(w)dudvdw = 
-oo J —oo J —oo 

/+oo r + CO r + CO 

p(u)du I p(v)dv I p(w)dw — \\p\\ 3 ■ 
-oo J—oo J — oo 



dx dudvdw = 



Corollary 3.1. It is clear that {Tp)(x) > for x G 3ft and Tp G L+(5ft). Also, it is 
straightforward to see that Tp is a continuous function, i.e. Tp G C(3?). 

Corollary 3.2. Consider the set B C LfQR) of PDF (Probability Density Func- 
tions), i.e. B = {p G : = 1}. It holds that if p G B then Tp G B. 
Clearly, if \\p\\ — 1 then \\Tp\\ = 1, i.e. the action of T on B conserves the number 
of particles in the gas. 

Corollary 3.3. Consider the following recursion: 

p n (x) = Tp n -i(x) with p G Lf($l) . 

Since \\Tp\ \ — \ \p\\ 3 therefore there are three possibilities depending on po: 

• If ||po|| < 1, then lini^oo ||p„(a,)|| = 0. 

• //| |po|| > 1, then lim^oo ||p„(x)|| = +oo. 

• V Iboll = 1, then \\p n (x)\\ = 1 Mn, i.e p n G B Vn. 
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Despite our interest resides specifically in the PDFs set, i.e. in B, some of the 
next theorems are also valid for generic elements of L\ . 

Theorem 3.2. The mean value of the velocity in the recursion p n = T n p is 
conserved in time. In fact, it is null for all n: 

< x, Tp >=< x, T 2 p >=< x, T 3 p >=•••=< x, T n p >=•••= , 

where 

P + OO 



/+oo 
f(x)g(x)dx . 
-oo 



It means that the zero total momentum of the gas is conserved in its time evolution 
under the action of T . 

Proof. Assume that p € B then 

p+oo 



/+oo 
xTp{x)dx = 
-oo 

' r ' p{u)p(v)p{w) 



SSL 



dudvdw I dx 



u 2 +v 2 +w 2 >x 2 2\Ju 2 + V 2 + W 2 J 

/+oo />+oo />+oo / \ / \ / \ / \/u 2 -\-v 2 +W 2 \ 
/ / rtwrfrQ ( / v xdx ) dudvdw = 
-oo J — oo J —oo 



" ' * ■• + °° p(u)p(v)p(w) ' 
y. J x J -oo 2\Ju 2 + V 2 + W 2 



2 -i + \/ u 2 +v 2 +w 2 

dudvdw = . 

-\/u 2 +v 2 -\-w 2 



Now, we know that if p € B then Tp € B, therefore < x, T n p >= Vn. □ 

Theorem 3.3. The total energy of the gas is conserved in time. In other words, 
for every p € B we have 

< x 2 ,p >=< x 2 ,Tp >=< x 2 ,T 2 p >=< x 2 ,T 3 p >= ■ ■ ■=< x 2 ,T n p >= ■ ■ ■. 
Proof. Assume that p € B then 

/+oo 
x 2 Tp(x)dx = 
-oo 

= r *>([[[ p^ph dudvdw \ dx= 

J-oo \J J Ju 2 +v 2 +w 2 >x 2 2Vu 2 + V 2 + W 2 J 

. r- p- r- pmpW)pm ( 1+*™*+ a iuAviw _ 

J —oo J —oo J —oo 2^/U 2 +V 2 +W 2 \J-^u 2 +v 2 +w 2 I 
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/+oo />+oo Z'+OO 
-oo J — oo J —oo 



p(u)p(v)p{w) 
2Vu 2 + v 2 + w 2 



dt/di;diu . 



+oo ,+oo ,+oo p^Jp^Jp^) 2 . 2 2 2 ,, , , 

- (u 2 + v 2 + w 2 ) 3 dudwdw = 



OO «/ — OO </ — OO 



2v^ 



oo />+oo />+oo 



(w 2 + w 2 + u> 2 )p('u)p(v)p(w;)d'ududw 



oo J — oo ^ — oo 



2 /-t-oo ^ /-t-oo 2 /-t-o° 

= — u 2 p(u)du + — v 2 p(v)dv + — / u> 2 p(u;)chz; = 

3 J— oo 3 J— co 3 J— oo 

/+oo 
x 2 p(a;)da; =< x 2 ,p > . 
-oo 

Now, wc know that if p g _B then 2> g B, therefore < x 2 , T n p >=< x 2 ,p > Vn. □ 

Theorem 3.4. TTie one-parametric family of PDFs p a (x) — y/^e~ ax2 , a > 0, are 
/jzed points of the operator T. In other words, Tp a — p a . 

Proof. When a = then p Q = 0, which is clearly a fixed point of T. Suppose now 
that a^O, 

Pa{u)p a (v)p a {w) 



Tp a 



2 +v 2 +w 2 >x 2 2\/u 2 + V 2 + w 2 



dudvdw = 



l l lu 2 



a rw e -a(u 2 +v 2 +w 2 ) 

7T V 7T 



dudwdw. 



a 2 +v 2 +w 2 >x 2 2\/u 2 + V 2 + W 2 

Transforming the integral region to spherical coordinates by the change of variables 
u = r smifjcos0,v = r sin?/> sin0, u; — rcostp, 



the proof becomes straightforward: 

hOO />7T r-2-K 



j—\-00 p7T i>27T j 

j\x\ Jo v ^ 



e ar rsmtpdOdipdr = 



■ e ar r sin ipdtpdr 



r+oo 

= 2a 

J\x\ 



re' ar dr = 





+co ! 






\x\ ~ V 



e- ax = Po 



□ 



Theorem 3.5. The operator T is Lipschitz continuous in B with Lipschitz constant 
S < 3. 
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Proof. Suppose p(x),q(x) G B, then 

\\T P - Tq \\ =r [ f f ^(y-^w dudvdw dx = 

J-oo J J Ju 2 +v 2 +w 2 >x 2 2 V U 2 + V 2 + W 2 

< ( +0 ° ( ( ( \p(u)p(v)p( W ) - q (u) q (v)q( W )\ dudvdwdx = 

J-oo J J J u 2+v 2 +w 2 >x 2 2\/lL 2 + V 2 + W 2 



oo •/ — oo •/ —CO 



Vi^+^+u^ 2\/w 2 + v 2 + u> 2 

/+oo Z-+00 Z> + 00 
-oo J —oo J —oo 

/+oo r+oo r+oo 
\ \ \p(u)p(v)p(w) - p(u)p(v)q(w) + p(u)p(v)q(w)- 
-oo J —oo J —oo 

—p(u)q(v)q(w) + p(u)q(v)q(w) — q(u)q(v)q(w)\ dududu> 

/+oo Z> + 00 r+oo 

p(u)du I p(v)dv / \p(w) — q(w)\dw+ 
-oo J —oo J —oo 

/+oo Z> + 00 Z> + 00 

p(u)du / q(w)dw / \p(v) — q(v)\dv+ 
-oo J —oo J —oo 

/+oo r+oo Z-+00 

q(v)dv / q(w)dw / \p(u) — q(u)\du = 
-oo J —oo J —oo 

/+oo 
\p(x) - q(x)\dx = 3\\p - q\\ . 
-oo 



□ 



Theorem 3.6. Suppose that lim„^co \\T n p{x) — n{x)\\ = 0, and is a contin- 
uous function in B, then p{x) is the fixed point of the operator T for the initial 
condition p(x) G B. 

Proof. First we prove that for each e > there exists M so that Vn > M: \ \T n p — 
T/i\ | < e. Since lim„_ ) . 00 \ \T n p(x)— fi(x)\\ = 0, we can say that 37V such that Vn > TV: 
\\T n ~ 1 p — p\\ < e/3. Now, we choose M = N . Then, due to the Lipschitz continuity 
of T in B, we have 

Vn > M : ||T>-7>|| < ^T^p - u|| < 3- | = e. 

Now, by uniqueness of the limit, it implies that Tp — /i. Therefore, for the initial 
condition p(x) G B, fi(x) is the fixed point of T. □ 
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Conjecture: For any p E B, with finite < x 2 ,p > and verifying 
liixin-yoo \ \T n p(x) — n(x)\ \ — 0, the limit fi(x) is the fixed point p a (x) = \f%^~ ax , 
with a — (2 < x 2 ,p >)~ . In physical terms, it means that for any initial velocity 
distribution of the gas, it decays to the Maxwellian distribution, which is just the 
fixed point of the dynamics. Recalling that < x 2 ,p >= kr, with k the Boltzmann 
constant and r the temperature of the gas, and introducing the mass m of the 
particles, let us observe that the MVD is recovered in its 3D format: 



MVD = 



3 

p a (u)p a (v) Pa (w) = (—J 2 e -™*(« 2 +« 2 +- 2 ) wi th a = (2k T y 



4. Examples 

In this section, we guess by simulation of many examples the truth of the above 
Conjecture, that is: 



lim 

n— f oo 



7T 



= with a = (2 < x 2 ,p >) 1 



T n P (x) - 
See the next examples. 

Example 4.1. Take p(x) a constant even function with mean value equal zero: 

I \ \\ ~ 1 < x < !) 
p(x) = < f 

[ otherwise, 
so a — (2 < x 2 ,p >) _1 = 3/2. Therefore we expect that 



T n p(x) — > fj,(x) = ^ — e- 
See the graphs of p(x), Tp(x) and /i(x) in Fig. [T] 




Fig. 1. Example l4.ll (a) p(x) and fJ.(x), (b) Tp(x) and p.(x). 
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Example 4.2. Take now p{x) a constant even function with mean value different 
from zero: 

p(*) = /* °<*< 2 > 
]_ otherwise, 

so a = (2 < x 2 ,p >) _1 = 3/8. Therefore we expect that 



T>(x) ^ /i(x) = y — e 
See the graphs of p(x), Tp{x) and /z(x) in Fig. [2J 




(a) (b) 
Fig. 2. Example 14.21 (a) p(x) and //(a;), (b) Tp(x) and //(a;). 



Example 4.3. Assume p(x) a non-constant even function with mean value equal 
zero: 

p{x) = ^iT^) ' 

so a = (2 <x 2 ,p>)~ 1 = l/2. Therefore we expect that 

r P (,) B -^ M (x) = £ e ^. 

V Z7T 

See the graphs off>(x), Tp(x) and /i(x) in Fig. [3l 

Example 4.4. Take p(x) a piecewise linear function with mean value equal zero: 

\ -1 < x < 0, 
p(x) = { \ < x < 2 , 
otherwise , 
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Fig. 3. Example 14,31 (a) p(x) and p-(x), (b) Tp(x) and p.(x). 



so a = (2 < x 2 ,p >) 1 = 3/5. Therefore we expect that 



See the graphs of p(x), Tp(x) and //(a;) in Fig. |4] 



o.s 



O.d 



0.4 



0.2 




Fig. 4. Example 14.41 (a) p(x) and /u(a:), (b) Tp(x) and 



Example 4.5. Finally, take p{x) a more complex shaped function with mean value 
equal zero: 



p(x) = —= x e 

\/7T 
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so a = (2 < x 2 ,p >) _1 = 1/3. Therefore we expect that 

T>(x) n -^° fa) = 
See the graphs of p(x), Tp(x) and ji(x) in Fig. [5j 




(a) (b) 
Fig. 5. Example 14.51 (a) p(x) and fi(x), (b) Tp(x) and fi(x). 



Let us observe that even if the initial II? distribution p(u) do not have the 
spherical symmetry in the 31? format p(u)p(v)p(w) (see for instance Example 14. 4[) . 
the operator T seems to recover this symmetry just after the first iteration. This 
suggests that in general the decay to the MVD is very fast, and after some few it- 
erations it would be difficult to distinguish how far is the system of the asymptotic 
Maxwellian regime (MVD). In fact, we have also checked the increasing of entropy 
with time (H-Theorem), that is verified in all these examples. This is another im- 
portant point now under study, that was also addressed in Boltzmann theorypQ and 
that can be the subject of a future work for the present model. 

5. Conclusions 

In this work, a new model for ideal gases has been presented. This model gives ac- 
count of the decay of the velocity distribution of an ideal gas toward the Maxwellian 
distribution. The first step has been to construct, under two simplistic and realistic 
assumptions, a discrete time evolution operator in the space of velocity distributions. 
The second step has consisted of showing the physical properties of this operator, 
specifically the conservation in time of the total energy and zero momentum of the 
gas. Third, the Maxwellian velocity distribution has been proved to be a fixed point 
of the dynamics toward which the gas asymptotically evolves with time. Finally, 
several examples, that display this behavior, have been explicitly given. 
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Let us conclude by saying that the type of models explained here have been 
successful in explaining the wealth distribution provoked by random markets and 
also the velocity distribution reached by ideal gases. We hope that more problems 
coming from other contexts can be studied under this new framework in the next 
future. 
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